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Beam  Trapping  in  High  Current  Cyclic  Accelerators 
with  Strong  Focusing  Fields 

I.  Introduction  y 

,«  2  3 

In  cyclic  induction  accelerators^  such  as  the  conventional  *  and  the 
modified  betatron^-^ -.the  energy  of  the  particles  increases  slowly  in 

synchronism  with  the  vertical  (betatron)  magnetic  field.,  For  example,  the  NRL 

\ 

O 

modified  betatron  has  been  designed  to  increase  the  energy  of  the  gyrating 

particles  by  -  1/3  KeV  per  revolution. 

As  a  consequence  of  the  slow  acceleration,  the  charged  particles  must  be 

confined  by  the  weak  focusing  magnetic  field  over  long  periods  of  time, 

Q  I  (V-  11  —  15 

typically  milliseconds,  and  thus  field  errors,  *  instabilities  and 

radiation  losses  can  impose  limitations  on  the  acceleration  process.  These 
limitations  can  be  substantially  relaxed  if  the  acceleration  were  to  occur 
rapidly,  say  over  a  few  microseconds.  An  appropriate  name  for  such  an 
accelerator  is  REBA-TRON  (Rapid  Electron  Beam  Accelerator). 

A  possible  configuration  for  a  rebatron  is  shown  schematically  in  Fig. 

1.  The  high  gradient  localized  electric  field,  responsible  for  the  rapid 
acceleration,  is  produced  by  a  convoluted  parallel  transmission  line,*'’ 
although,  other  transmission  lines  are  also  suitable.  Since  the  acceleration 
occurs  over  a  few  microseconds,  the  vertical  magnetic  field  cannot  be 
increased  in  synchronism  with  the  energy  of  the  particles.  Therefore,  a 
strong  focusing  field  is  needed  to  confine  the  high  current  electron  beam. ^ 

The  confining  properties  of  torsatron  magnetic  fields  in  rebatrons  have  been 

21 

studied  recently  and  the  results  are  very  promising. 

In  this  paper  we  consider  a  possible  mechanism  which  could  trap  a  high 

current  electron  beam  in  the  strong  focusing  magnetic  fields  of  the  rebatron. 
Manuscript  approved  January  18,  1985. 
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The  beam  trapping  process  in  rebatrons  is  inherently  difficult  since  the  strong 

focusing  fields  result  in  particle  orbits  which  are  relatively  insensitive  to 

20 

an  energy  mismatcn.  Thus  small  changes  in  the  betatron  magnetic  field  are 
not  sufficient  to  move  the  beam  from  the  injection  point  to  the  center  of  the 
chamber. 

The  proposed  injection  scheme  is  based  on  the  presence  of  a  drag  force 

that  acts  on  the  beam  centroid  when  the  chamber  walls  are  resistive. AJ 

Consider  a  pencil-like  electron  beam  inside  a  straight,  finite  conductivity 

chamber.  In  the  absence  of  strong  focusing  fields  the  image  or  induced 

forces,  due  to  the  self  fields,  cause  the  beam  centroid  to  undergo  transverse 

motion  within  the  chamber.  In  a  perfectly  conducting  chamber  the  induced 

force,  F,  which  is  directed  normal  to  the  chamber  wall,  together  with  the 

external  longitudinal  magnetic  field  produce  a  F  x  B  force  on  the  beam.  This 

force  causes  the  beam  to  execute  transverse  bounce  oscillations  within  the 

chamber.  However,  when  the  resistivity  of  the  wall  is  finite,  the  decay  of 

the  wall  currents  produce  a  component  of  induced  force  which  is  directed 

tangentially  to  the  wall.  For  a  straight  electron  beam,  this  drag  force, 

together  with  the  longitudinal  magnetic  field,  causes  the  beam  to  spiral 

outward  towards  the  chamber  walls.  For  a  cyclic  electron  beam,  the  beam 

2  2 

centroid  spirals  inward  towards  the  center  of  the  chamber  if  n  r,  /a  <  1/2, 

s  b 

where  ng  is  the  self  field  index,  r^  is  the  beam  radius,  a  is  the  wall  radius 
and  the  external  field  index  has  been  set  equal  to  1/2. 

In  the  presence  of  strong  focusing  fields,  the  bounce  (poloidal)  motion 
of  the  beam  is  due  mainly  to  the  gradients  of  the  external  magnetic  fields  and 
its  direction  can  be  selected  to  provide  an  inward  drift  velocity  of  the  beam 
centroid,  resulting  in  an  inward  spiraling  motion. 


The  implications  of  the  resistive  wall  on  the  stability  of  the  beam  must 
be  carefully  considered.  However,  since  the  relativistic  cyclotron  frequency 
corresponding  to  the  vertical  field  varies  rapidly  with  time,  it  is  likely 
that  longitudinal  instabilities  and  in  particular  the  resistive  wall  mode  can 
be  avoided. 

The  present  injection  scheme  is  probably  also  applicable  to  the  Racetrack 

o  9 

Induction  Accelerator,"'"  a  device  similar  to  the  rebatron  but  without  a 
vertical  magnetic  field. 

To  determine  how  rapidly  the  beam  can  be  moved  from  the  wall  to  the  axis 
of  the  chamber,  we  have  analyzed  the  beam  dynamics  inside  a  straight 
cylindrical  pipe  having  resistive  walls.  Our  analysis  indicates  that  for  a 
fairly  wide  range  of  parameters  the  beam  centroid  can  spiral  inward  from  its 
injection  position,  near  the  wall  of  the  chamber,  to  the  axis  of  the  chamber 
in  less  chan  50  nanoseconds. 


II.  Dynamics  of  the  Beam  Centroid 

To  make  the  problems  analytically  tractable,  the  dynamics  of  the  beam  is 
studied  in  cylindrical  geometry,  i.e.,  toroidal  corrections  associated  with 
the  fields  are  neglected. 

Consider  a  pencil-like  electron  beam  inside  a  straight  cylindrical  pipe 
of  circular  cross  section  as  shown  in  Fig.  2.  The  pipe  thickness  is  denoted 
by  A,  the  minor  radius  by  a,  and  the  wall  conductivity  by  a.  The  centroid  of 
the  beam  is  located  at  the  coordinates  (x,y)  relative  to  the  center  of  the 
minor  cross  section  of  the  chamber. 

The  transverse  components  of  the  torsatron  fields  or  rotating  quadrupole 
magnetic  field  near  the  z-axis  are  of  the  form 
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(la) 


b  =  b  k  (y  cos  k  z  -  x  sin  k  z), 
X  0  w  w  w 


b  =  b  k  (x  cos  kz+vsinkz),  (lb) 

vow  w  w 


where  k^  =  2tt/jI,  £  is  Che  period  of  the  magnetic  field  and  bQ  is  a  constant 
that  is  a  measure  of  the  field  strength.  Along  with  the  external  transverse 
periodic  field  in  Eq.  (1)  we  also  have  a  constant  uniform  axial  magnetic  field 
given  by  B  =  B  e  ,  which  is  the  sum  of  the  torsatron  axial  field  and  any 
additional  external  axial  magnetic  field.  The  induced  electric  field  due  to 
image  charges  on  the  walls  of  the  conducting  chamber  is  not  affected  by  the 
finite  wall  conductivity  and  is  given  by 


where  a  is  the  radius  of  the  chamber,  mQ  is  the  electron  mass,  c  is  the 

velocity  of  light,  |e|  is  the  magnitude  of  the  electronic  charge,  and  v  is 

2 

Budker's  parameter.  This  parameter  can  be  expressed  as  v  =  (cu,r,/2c)  , 

b  b 

2  1/2 

where  *  (4ir|e|  mo/mQ)  is  the  beam  plasma  frequency  and  rb  is  the  beam 
radius.  The  induced  magnetic  field  at  the  center  of  the  beam  is  due  to  image 
currents  flowing  in  the  finite  conducting  wall.  This  magnetic  field  is 
modified  by  the  finite  wall  conductivity  and  can  be  written  in  the  form 


~ind 


2m  c 
o 


— 2"  SQ  G(t-x)  f  y(x )  ex  -  x(t)  ey)  dr. 


(3) 


where  3^  *  v0^c»  vo  c^e  constant  axial  beam  velocity  and  G(t-x)  is  the 
magnetic  diffusion  function.  The  induced  magnetic  field  in  Eq.  (3),  in  the 
form  of  a  convolution  integral,  denotes  the  fact  that  the  field  is  a  nonlocal 


4 


function  of  time.  The  function  G(t — has  been  derived  elsewhere^  and  we 

will  simply  state  it  when  needed.  For  infinitely  conducting  walls 

u(:-t)  *  o(c-t)  and  the  induced  magnetic  field  in  Eq.  (3)  reduces  to  its  usual 

form. 

The  transverse  velocity  of  the  beam  centroid  is  nonrelativistic,  i.e. , 

•  * 

i x | ,  |y|  «  c,  whereas  the  axial  beam  particle  velocity  is  highlv  relativistic 
and  is  assumed  constant.  Using  the  external  field  in  Eq.  (1)  and  the  induced 
fields  given  in  Eqs.  (2)  and  (3),  together  with  the  external  solenoidal  field 
®oez’  t'le  equations  governing  the  transverse  motion  of  the  beam  centroid  are 


x(t)  -  u)  x(t)  +  Q  y(t)  -  to  a)  !x(t)  cos  tu  t  +  y(t)  sin  m  t| 
s  o  o  w  w  WJ 


III  t 

1  <5G(  t-r )  x(t)  dx, 

3  0  0  0 


y(t)  -  a)  y(t)  -  Q  x(t)  +  a)  to  I  y ( t )  cos  u  t  -  x(t)  sin  w  t' 

Q  tJ  1  J  r.r  r.T  J 


III  t 

w  1  5  G(t-x)  y(x)  dx, 

s  o  o  o 


l  l  i  l 

where  u>  =  2c  v/(y  a  ),  ft  =  |e|  B  /(y  m  c),  u  =  Ie|b  /(y  m  c),  u  =*  k  v  , 
s  o  o  ooo  o  ooo’wwo 

,  2  "1/2 

yQ  =  i  1“8  I  and  5G(t-x)  =  G(t-x)  -  S(t-x).  The  right  hand  side  of  Eqs. 


(4a, b)  contains  the  finite  wall  conductivity  effects  and  vanishes  for 
perfectly  conducting  walls.  In  obtaining  Eqs.  (4a, b)  we  have  set  z  =*  v^t  in 
Eqs.  (la,b)  since  we  are  seeking  the  time  evolution  of  a  cross-section  of  the 
beam  centroid.  Rather  than  solving  Eqs.  (4a, b)  directly  it  has  been  found 
convenient  to  first  transform  these  equations  to  a  frame  rotating  with  a 
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F  (5  ,  t*o)  'j(s), 


Fv(s,t=o)/D(s) , 


where  x(s),  v(s)  are  the  transforms  of  x(t)  and  y(t)  respectively, 


2  2  2  2  2  *.  2  2  222*. 

D(s)  =  {  s  +  Q  +  u  y  S  5C(  s)  1  f  s  +  ft  +  u  f  B  5  C  (  s)  ' 
xsco  vsoo 


2  2  2  *.  2 

+  '  (ft  -(jj)c+a)Yg5S(s)) 

o  w  s  o  o 


5C(s)  =  7  (  6G( s  -  iw  /2)  +  «$G  (s  +  ioj  /2)!, 


5S(s)  =  77  ’  $G(s  -  i»  /2)  -  <5G(  s  +  ioj  /2 ) ) 

4-1  W  W 


are  the  transforms  of  6G(t)cos(<D  t/2)and  oG(t)sin(oj  t/2)  respectively. 

w  w 

The  initial  conditions  on  x,  y,  x  and  v  are  contained  within  the 
functions  Fx(s,t=o)  and  Fy(s,t=o).  The  exact  forms  of  these  functions  are  not 
of  interest  since  we  will  be  mainly  concerned  with  the  roots  of  the  function 
D(s)  which  is  independent  of  initial  conditions. 

For  an  infinitely  conducting  chamber,  i.e. ,  when  o  *  « ,  the  functions 
5C(s)  =  5S(s)  vanish  and  D(s)|  =  D  (s)  becomes 


2  2  2  2  2  2 

D(s)=(s+ft)(s  +  ft  )  +  s  (ft  -  u>  )  , 

*  x  y  o  w 


the  roots  of  D  (s)  *  0,  which  will  be  denoted  bv  s.*,  determine  the 

oo  J  O 

characteristic  frequencies  associated  with  the  beam  centroid  motion  and  are 


given  by 
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.i  +  ( r.  -  i)  )  1  , 

-Z-. - 2 - *  i  ;  (a  +  ,1 


2 

2 

x 


O'  - 

o 


2 

w  )  > 


2  2 
"x  v 


1  '? 


(ID 


Performing  a  Taylor  expansion  about  s  -  sQ  we  find  that  the  roots  of 
D(s)  =  0  in  Eq.  (9)  are  given  by 


D(s  ) 

_ o 

S  ~  s0  3D(s)/Ss| 


(12) 


s  =  s 


where  s  is  given  in  Eq.  (11), 


D(s  ) 
o 


2  2  2  2 
2'jj  v  8  !  i  s 

s  o  o  l  o 


2  2 

n  +  q 

x  y 

o 


«G(s  ) 
o 


(n  -  ta  )  s  «  S(s  )!  , 

O  WO  0  1 


(13a) 


and 


3D(s) 

3s 


s=s 


2  2  2  2 

4s  f  s  +  (S3  +  S3  +  (ft  -  u  )  )/2 ! . 

oo  x  y  o  w 


(13b) 


Substituting  Eqs.  (13a, b)  together  with  Eqs.  (9b, c)  into  Eq.  (12)  yields 


2  2  2.  2  2 
w  y  3  1 2,  ( s  )6G(s  +  iu  12)  +  S3  (s  )<SG(s  -  lu/2) 

s  o  o  1  +  o  o  w  -o  0  WJ 

3  =  So - 2 - 2 - 2  2 


(14a) 


4s  [  s  +  03  +  S3  +  (n  -  to  )  ) / 2 ' 

o  ‘  o  x  v  o  w  J 


where 


2  2  2  2 

H  ( s  )  =  s  +03  +  Q  )/2  +  is(fl  -  «). 

4-0  o  x  v  -  o  o  w 


(14b) 


As  mentioned  earlier,  the  derivation  of  the  Laplace  transform  of  the 

diffusion  term,  ■SG(s)  =  G(s)-1  will  not  be  repeated  here.  It  has  been  derived 
1  1 

elsewhere  and  found  to  be  given  by 
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5G(  s) 


2/(F(s)-l), 


( 15a) 


Z 

5 

r 

t- 

k 

I 

! 


where 


[ 1  +  7  u(s)  tanh  y  (s)  j 

F(s)  =  -  —  u  ( s )  - r - 

[tanh  y(s)  +  ^  y  (s)  ] 


and 


y(s)  =  A  [4iras/c  + 


2  1/2 

1/a  ) 


(15b) 


(15c) 


« 

m 

I 

» 

. 

► 


Equation  (14)  together  with  Eq.  (15)  completes  the  formal  derivation  for  the 
roots  of  D(s)=o  which  in  turn  determines  the  dynamics  of  the  beam  centroid 
within  the  resistive  chamber. 

The  real  part  of  s  determines  the  macroscopic  motion  of  the  beam.  If 
Re(s)  <  0  the  beam  centroid  spirals  inward  towards  the  center  of  the 
chamber.  If  Re(s)  >  0  the  beam  is  unstable  and  its  centroid  spirals  outward, 
away  from  the  center  of  the  chamber.  To  utilise  this  mechanism  for  trapping 
the  beam  it  is  required  that  Re(s)  <  0  and  that  the  magnitude  of  Re(s)  be 
sufficiently  large  in  order  for  the  beam  to  spiral  to  the  center  of  the 
chamber  in  a  reasonably  short  period  of  time. 


Ilia.  Growth  Rate  in  the  Absence  of  Torsatron  Fields 


We  first  consider  the  trivial  case  where  the  torsatron  magnetic  field 


vanishes,  i.e.,  to  *  0.  In  the  absence  of  the  strong  focusing  field  it  is 


unnecessary  to  transform  to  the  rotating  frame,  hence,  we  also  take  to  *  0. 
2  2 

For  to  /fj  <<  1  we  find  that 
s  o 


9 


3  =  3  + 

O 


a)  y  3  5G(s  ) 
s  o  o  o 


s  2 
o  o 


(16) 


wh 


ere  =  ito^  and  uig  =  ±  cns/  f 1  *  The  real  Part  of  s  is  given 


bv 


h  2  2 
U3  Y  8 

Re(s)  =  T  =  —  0/  £m6G(iajB)  , 


Vo 


(17) 


where  <5G  is  given  by  Eq.  (15). 

When  the  chamber  thickness  A  is  much  greater  than  the  skin  depth 
associated  with  the  beam  centroid  bounce  frequency,  i.e.  ,  »  1,  we  find 

O 

that 


<5G?  s  =  iuu  =  +  irn  /2  ) 
v  o  B  ~  so 


-  2^2  ( cos  w/4  -  i  -j^j— r  sin  tt/4)  ,  (18) 


where  <5fi  =  c/(2tt<jI  I  ) 1 1 2  is  the  skin  depth.  Substituting  Eq.  (18)  into  (17) 
yields 
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which  always  results  in  instability. 


Illb.  Decay  Rate  in  the  Presence  of  Torsatron  Fields 

In  the  presence  of  the  strong  focusing  field,  i.e.,  when  u>  *  0,  we  may 

o 

L 

neglect  the  self  field  contribution,  i.e.,  o»g ,  in  the  evaluation  of  the 
characteristic  frequencies  of  the  macroscopic  beam  motion. 


Using  the  definitions 
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and  assuming  ew,  eQ  <<  1  we  find  that  the  characteristic  frequencies  given  by 
Eq.  (11)  are 
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(21) 


Since  «  1,  the  upper  sign  in  Eq.  (11)  represents  the  slow  drift  motion 

of  the  beam  centroid  and  the  lower  sign  denotes  the  fast  motion.  In  what 

follows  only  the  slow  macroscopic  beam  motion  will  be  considered  since  we  are 

primarily  interested  in  beam  trapping. 

From  Eq.  (21),  the  roots  of  D(s)|  “0  corresponding  to  the  slow  beam 

o  a  « 

motion  are  given  by  sq  *  iuig,  where 
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(22) 


2 

is  the  bounce  frequency.  Since  <<  1,  the  bounce  frequency  in  the  rotating 
frame  is  approximately  (d^/2. 

For  >  0  the  roots  of  D(s)  *  0,  given  in  Eq.  (14a),  become 
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To  evaluate  the  real  oart  of  s  in  Eq.  (23)  we  require  the  imaginary  part 
of  <$G( i( uu  +  iUg/2)),  where  5G  is  given  by  Eq.  (15).  Using  Eq.  (15c)  it  is 


found  that 


(  (  210-0^/4)  )1/2 
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where 


6  =  c  (2iro|  oj^I  j 


is  the  skin  depth  associated  with  the  frequency  u)  *  v  k  . 
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Case  ( i) ,  A/ 6  «  1 

For  A/6  «  1,  <5G  is  given  by 


60(1(0^  *  o>w/2)) 


and  the  decay  rate  becomes 
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where 
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'  =  aA/o  ,  e  =  w  Q  ,  e„  =  w_/H  ,  £,£„<<!  and  a/6  «  1. 
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Case  (ii)  ( A / 5 ) e  //2  >>  1 
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For  this  case  6G  is  given  by 
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and  the  decay  rate  becomes 
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In  both  (28)  and  (30)  we  have  made  the  assumption  that  e  ,  e  .  e  , 

s  w’  o' 

2 

e  <<  1  and  y  »  1. 

D  O 

The  decay  rates  given  by  (14)  together  with  (15)  are  shown  in  Figs.  (3) 
and  (4)  for  various  sets  of  parameters. 


Discussion 

In  this  paper  we  have  investigated  a  possible  mechanism  for  trapping  an 

intense  relativistic  electron  beam  confined  by  strong  focusing  fields.  In  our 

model  the  electron  beam  is  assumed  to  be  injected  into  the  torsatron  fields 

off  axis,  near  the  chamber  walls.  The  finite  resistivity  of  the  walls  results 

in  a  drag  force  on  the  beam  centroid  which  may  cause  the  beam  to  spiral  inward 

towards  the  axis  of  the  chamber.  We  have  analyzed  this  mechanism  and  obtained 

decay  rates  for  the  inward  spiraling  motion. 

As  an  illustration  of  the  trapping  process,  we  consider  a  15  kA(v  *  0.86) 

2  MeV(y  -  5)  electron  beam  injected  near  the  wall  of  a  10  cm  radius  resistive 
o 

chamber  of  thickness  A  ■  a/150.  The  axial  magnetic  field  is  taken  to  be 
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3  =  2.5  kG,  the  torsatron  coefficient  bQ  is  250  G  and  the  magnetic  oeriod 

l  of  the  torsatron  field  is  12.8  cm.  With  these  parameters  we  find  that 

e  =  0.02,  e  =  0.25,  e  =  0.2,  e„  =  0.12  and  0  =  8.8  x  109sec_1.  The 

s  w  o  B  o 

maximum  decay  rate  for  these  parameters  occurs  (See  Fig.  4)  at  A/8  =0.5  and 

-3  7-1 

is  T  *  5  x  10  SI  =  4.4  x  10  sec  .  The  conductivity  of  the  chamber  wall 
o 

oo  12-1  9-1 

is  a  ■  c^/(2ir|o)  |6“)  =  4  x  10  sec  ,  u>  =  e  0  =  2.2  x  10  sec  and 

w  w  w  o 

5  =  2A  =  0. 13  era.  The  e-fold  time  for  the  inward  spiraling  motion  of  the  beam 

is  t  =  1/F  -  22  nsec,  which  is  sufficiently  fast  to  trap  the  beam.  It  should 
e 

be  noted  that  in  this  example  for  l  =  12.8  cm  the  linearized  fields  of  Eq. 

w 

(1)  are  valid  only  within  a  small  distance  from  the  axis.  In  practice,  the 
more  exact  expressions  for  the  fields  should  be  used  and  the  problem  addressed 
numerically. 
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Schematic  of  a  rebatron  accelerator 


Conductivity  <r 

Figure  2  Model  of  cross  section  of  rebatron  chamber  showing  the  tranverse 

location  of  the  electron  beam.  The  beam  transverse  dynamics 
is  influenced  by  the  torsatron  field,  axial  magnetic  field  and 
induced  fields. 
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Normalized  Decay 


0.5  1.0  1.5  2.0  2.5  3.0 

A/8  (Thickness /Skin  Depth) 


Figure  3  Shows  the  decay  rate  for  the  beam  centroid  as  a  function 


Normalized  Decay  Rate,  T /  £1 


of  A/ <5  for  eq  =  0.1,  0.3,  0.6,  a/ A  =  100,  *  0.25, 
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